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Abstract. - We consider a massive inelastic piston, whose opposite faces have different coefficients 
of restitution, moving under the action of an infinitely dilute gas of hard disks maintained at a fixed 
temperature. The dynamics of the piston is Markovian and obeys a continuous Master Equation: 
however, the asymmetry of restitution coefficients induces a violation of detailed balance and a 
net drift of the piston, as in a Brownian ratchet. Numerical investigations of such non-equilibrium 
stationary state show that the velocity fluctuations of the piston are symmetric around the mean 
value only in the limit of large piston mass, while they are strongly asymmetric in the opposite 
limit. Only taking into account such an asymmetry, i.e. including a third parameter in addition 
to the mean and the variance of the velocity distribution, it is possible to obtain a satisfactory 
analytical prediction for the ratchet drift velocity. 
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Introduction. — Granular materials have been the 
subject of intense research in the last 20 years in 
Physics [1]. Most of the non-trivial phenomena that can 
be observed in a shaken box of sand are due to the in- 
elasticity of collisions among grains [2]. Kinetic energy 
is dissipated into heat, introducing an intrinsic time irre- 
versibility in the "microscopic" dynamics which can have 
consequences at a more macroscopic level: for instance 
species segregation [3], breakdown of energy equiparti- 
tion [4], apparent Maxwell-demon-like properties such as 
heat currents against a temperature gradient [5] or mass 
current against a density gradient [6] , ratchet-like net drift 
of a asymmetrically shaped tracers [7] [8] , rectification of 
thermal fluctuations [9] [10], and so on. It is well known, 
moreover, that the space asymmetry of an adiabatic pis- 
ton results in a stationary macroscopic motion [11]. Here 
we consider, instead, a model of asymmetric granular pis- 
ton with different coefficients of restitution which, with 
respect to a previously presented model of granular Brow- 
nian ratchet [7], has the advantages of being simpler to 
be studied analytically as well as realized in the labora- 
tory, and at the same time displays unexpected peculiar 
properties: in particular we will show how its stationary 
state is characterized by asymmetric velocity fluctuations, 
and that this asymmetry becomes crucial when the piston 



mass is smaller than the mass of surrounding disks. Such 
a light piston limit, which could be regarded as purely aca- 
demic in the case of a sect separating two molecular gases, 
is instead realistic in the case of granular gases, where the 
role "molecules" is played by large grains. 

Model. — The 2D model consists of a piston of mass 
M and height L surrounded by a dilute gas of N hard disks 
of mass m and density p. The faces of the piston have two 
different values of inelasticity, characterized by the coeffi- 
cients of restitution a±, on the left, and ai, on the right. 
The piston can only slide, without rotating, along the di- 
rection x, perpendicular to its faces. The thickness of the 
piston can be neglected for the purpose of our simulations 
and calculations. The particles-piston binary collisions are 
described by the rule: 



V 



V'-(l + a i )-—^(V'-v' x ) 
1 + 



v' x + (l + a i ) I ^(V'-v' x ) 



(1) 



where v and v' are the post-collisional and pre-collisional 
disk velocities respectively, while V and V are the corre- 
sponding velocities of the piston, while e 2 = m/M . Be- 
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cause of the constraint on the piston, its vertical velocity 
is always 0. The energy of this granular system is not 
conserved and an external driving mechanism is needed 
to attain a stationary state. In our model, the surround- 
ing gas is coupled to a thermal bath which keeps constant 
its temperature: the exact nature of the thermostat is 
not discussed here (many models have been introduced in 
the literature, see for example [12,13]), since we focus on 
the dynamics of the piston which is assumed not to cou- 
ple directly with the thermostat, but only with the gas 
particles. For the sake of simplicity we assume that the 
surrounding gas is homogeneous in space and time, with 
a Maxwell-Boltzmann velocity probability density func- 
tion (pdf) (f)(v x ,Vy) with zero average and given variance 
( v x) = ( v y) ( see below). Such an assumption can be con- 
sidered realistic as long as the gas is dilute and the char- 
acteristic time of coupling with the external heat bath 
is much shorter than that with the piston. In this case 
also Molecular Chaos for piston-disks collisions can safely 
be assumed, allowing for the use of the Direct Simulation 
Monte-Carlo (DSMC) algorithm to simulate the piston dy- 
namics [14]. Moreover, without loss of generality, we can 
choose one side of the ratchet to be elastic. Based on the 
rules (1), in fact, a system with mass M and coefficients 
of restitution < a 2 < ct\ < 1 is equivalent to one with 
ai = 1 an d effective parameters M' and < a' 2 < 1 given 
by the following relations 



M' 



m(l - cvi) + 2M 



1 



1 + ai 

( Q!l - Q.1 
1 + Q'l 



(2) 



In the following, therefore, we choose a.\ — 1 and we study 
the behavior of the piston for different values of the mass 
ratio e 2 and coefficient of restitution a 2 . 

Results and discussion. — Our choice of the initial 
conditions for the piston are V(0) = and X(0) = After 
an initial transient whose duration depends on all con- 
trol parameters such as collision frequency, coefficients of 
restitution, masses, etc., the piston reaches a stationary 
regime. In the following we shall indicate as T g = m(v^) 
the value of the gas temperature and T r = M((V— (V^) 2 ) 
the ratchet temperature. The DSMC simulations have 
been performed using T g = 1, m = 1. We focus here 
on the behavior of the piston, whose position and veloc- 
ity at time t are denoted as X(t) and V(t) respectively. 
Trajectories of the piston in its non-equilibrium stationary 
state, i.e. discarding transients, for particular choices of 
the parameters M and a 2 , are displayed in Fig. la. When 
the system is totally elastic or ol\ = a 2 no average mo- 
tion occurs. On the contrary, if the inelasticity of the left 
and right side are different, the piston shows a mechanism 
of rectification of disorder-induced fluctuations driven by 
the collisions and it drifts with average velocity (V) ^ 0. 
In particular this drift appears to be always oriented to- 
wards the side of the piston with smaller restitution coef- 
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Fig. 1: Panel a: Single trajectories, in DSMC, of the massive 
piston with e 2 = 0.01 (full symbol) and light piston with e 2 = 
10 (open symbol) for two different coefficients of restitution 
ai = 0.9 (circles) and ai = 0.3 (squares). Two cases with 
same inelasticity (ai — ai — 0.3) and different mass (e 2 = 
0.01 and e 2 = 10) are shown (diamonds). The trajectories in 
these cases are averaged over 100 realizations. Panel b and c: 
The rescaled probability density functions, in DSMC, for the 
different mass of the piston: e 2 = 0.01 (circles) e 2 = 1 (squares) 
and e 2 = 10 (diamonds) for ai = 0.3 (panel b) and ai = 0.9 
(panel c). The velocity V of the ratchet is rescaled by thermal 
velocity Vth = \jT r jM. The full line in the panel c represents 
a gaussian fit of the data. 



ficient, the positive direction in our case, so that (V) > 0. 
This phenomenon can be understood by considering the 
average momentum transferred by the gas to the piston, 
M(V' — V). Assuming that the piston is slower than the 
thermal velocity of the gas it is easy to show that it will 
experience a viscous drag force —jV, where 7 is given by 
(2 + ai + a 2 )py / M/(M + m)y/2mT g /'K plus a net force 



p M 
p — L-T 

~ 2 9 M + ; 



-(ai - a 2 ) 



(3) 



In other words, the momentum transferred by the gas is 
smaller on the more inelastic side, originating the drift 
in that direction. From the simulation one also observes 
that the velocity of the drift is bigger if the elasticity is 
smaller, while it seems not much influenced by the piston 
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mass. Mass, on the contrary, influences fluctuations in the 
trajectories, which are larger the lighter is the ratchet. In 
order to characterize the behavior of these fluctuations, we 
show, in Fig. lb and Fig. lc, the pdf P(V) of the rescalcd 
velocity of the piston. When the piston is massive e 2 <C 1, 
the distribution P(V) does not display any appreciable 
difference with respect to a Gaussian with finite mean. 
On the contrary, when the mass M of the piston is equal 
or smaller than that of surrounding disks, P(V) results 
asymmetric with a larger tail for positive velocities. Such 
an effect occurs at smaller elasticity and it is stronger if 
e 2 3> 1. Note also that in all our simulations we always 
have (V) > 0, while in some cases the maximum of the 
velocity pdf occurs at negative values. 

In the dilute gas limit, it is possible to study the piston 
dynamics by means of a Master Equation (ME) for P(V, t) 
which can be written as [7]: 



dP(y, t) 

dt 



dV [W{V\V')P(V',t) - W(V'\V)P(V,t)] 

(4) 



where the transition rate is: 

W(V\V) = pL^-lf [ dvQ[(-l) k (V - v x )] 
fc=i J 

(V - v x )</>(v)S[V — V' + fx k (V' - v x )] (5) 

with /i fe = (1 + a k )e 2 /(l + e 2 ) (see Eq. (1)) and 6 is the 
Heaviside step function. 

In order to obtain quantitative informations about the 
velocity and granular temperature of the ratchet we can 
write a system of equations describing the evolution of the 
first moments of the distribution, starting from (4) and 
(5), and invoke some approximation to obtain a closed 
set. By direct experience we have found that, due to the 
asymmetry of P(V), it is essential to consider a further 
parameter £ = ((V — (V)) 3 ) in addition to the average 
velocity and granular temperature [15]. To this aim we 
assume that the pdf of the piston can be written as 



P(V) 



71=0 



exp 



O) 



Imposing the normalization and that ((V — (V)) 2 ) = 
T r /M, the series (6), with a truncation at n = 3, becomes 



retaining only the terms of first order in (V), one obtains 
the following differential equations 



P M 



3 V 27rT r 



03(7?)? 



p / M 3 



6 V 2nT r 



(8) 



(9) 



dd 3T 2 p / 2T3 



dt 



2 V irM 3 



(10) 



where r] = T r /T g . The coefficients of the above equations 
arc given in the appendix. 




Fig. 2: The velocity of the piston, rescaled with its thermal 
velocity Vth = \fT r /M, as function of the coefficient of the 
restitution 012 for different values of the parameter e 2 : 10 (cir- 
cles), 1 (squares), 0.1 (diamonds) and 0.01 (triangles). The 



P(V) 



M 
2vrT r 
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6dV 3 



exp 



M , 



In the above expression £ represents the first term that is 
a measure of the asymmetry of P(V) about the average 
value. Under these assumptions the equations for (V), T r 
and £ can be obtained multiplying both sides of Eq.(4) by 
V, M(V- (V)) 2 and (V- (V)) 3 respectively and perform- 
ing the integrations. After some calculations, expressions 
in analytical form for these equations can be derived if it 
is assumed that (V) <C Vth, where Vth = \jT r jM. By 



symbols correspond to the simulation data while the lines cor- 
respond to the solutions obtained from Eqs.(ll)-(13). Inset: 
(V)) 2 (7) The velocity of the piston, rescaled with the thermal velocity 
of the gas v g th = \/T g /m, as function of the coefficient of the 
restitution Q2 for the same values of e 2 . The dashed and dot- 
dashed lines correspond to the theoretical calculation for the 



cases e = 10 and 6=1 respectively, assuming P(V) exactly 
gaussian. 



We can solve Eqs.(8)-(10) in the stationary state (all 
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Fig. 3: The temperature ratio rj — T r /T g as function of the 
coefficient of the restitution a? for the same cases of Fig. 2. 
The symbols correspond to the simulation data while the lines 
correspond to the solutions obtained from Eqs.(ll)-(13). 




Fig. 4: The third moment £ = {(V - (V)) 3 }, rescaled with the 
quantity Vg th as function of the parameter e 2 for a?2 = 0.3. 
The symbols correspond to the simulation data while the line 
corresponds to the solution obtained from Eqs.(ll)-(13) 



time derivatives are put to zero), obtaining: 

2 V 2 * 302(77)534(77) + 03(77)324(77) 
3T/ 2 

203(77) 

3V2^V th [2b 2 ( v )V t l(V) + 26 3 C] + 2b i (r l )(V}£+ 

24^(77) = (13) 

where 314(77) = 64(77)01(77) + 461(77)04, 324(77) = 
64(77)22(77) + 3^62(77)04 and 334(77) = 64(77)03(77) - 67r6 3 c 4 
(see appendix). 

A numerical solution of system (11)-(13) gives the sta- 
tionary values of (V), 77 and £. In the symmetric case 
a i — a 2 — Q!, it turns out that 01(77) = 62(77) = 63 = 
01(77) = c 4 = the system has the solutions (V) = £ = 
and from (13) we get 77 = (1 + a)/[e 2 (l — a) + 2], that is 
the same solution already obtained by Martin et al. [16]. 
When ai = oli = 1 the temperature of the piston is equal 
to that of the gas. The results in the asymmetric case are 
shown in Figs.(2)-(4) together with the simulation data. 

Fig. 2 displays a good agreement between theory and 
DSMC for the rescaled observable (V) /Vth for all values 
of «2 and e 2 studied. It is remarkable that our results are 
good even if the above assumption does not hold true, i.e. 
when (V) ~ Vth- In the inset of Fig. 2 we show the behav- 
ior of (V) I Vgth as a function of the coefficient of restitution 
02 for different e 2 . The theoretical results for [M < m), 
derived assuming a Gaussian P(V) (i.e. £ = 0) are also 
displayed. They appear to disagree with the numerical 
results, suggesting that in these cases a Gaussian pdf for 
the piston velocity fluctuations is a very poor assumption. 
The third moment £, that is a marker of the asymmetry of 
the pdf, is not small (see Fig. 4) and expansion (6) is cru- 
cial. The small differences between theory and DSMC for 



e 2 = 10 (more visible in the inset of Fig. 2) are probably 
due to the fact that we have considered only the first three 
terms of the expansion. In fact the numerical results show 
for large values of e 2 the fourth moment becomes relevant 
so that one should take into account the evolution of such 
a quantity in the coupled set of equations in order to give 
more accurate results. 

We would like to add a comment concerning the high 
velocity tails of the distribution function P(V). In spite 
of the fact that the body of such a distribution is not 
Gaussian and skewed, we have been able to prove that its 
tails are of Gaussian nature, but with different effective 
temperatures in each tail, reflecting the asymmetry of the 
piston. In particular, the temperature of each tail is given 
by a formula similar to the one presented by Martin and 
Piasecki for the inelastic intruder: 



Ti = T 



1 + a, 



(1 - a t )m/M + 2 



(14) 



Finally, we have studied, in Fig. 5, the diffusion coef- 
ficient D of the piston for different values of e 2 . In order 
to better characterize the situation of a light ratchet, we 
have calculated D in two different ways. First we have 
measured it by its definition D = limt_ ) . 00 IMllz^MlilJ ) 
where d(t) = X(t) — X(Q) is the displacement of the tracer 
with respect to its initial position at a time t = 0, taken 
when the whole system has become stationary. Second we 
have used the self-correlation function of the piston veloc- 
ity defined by C(t) = ((V(t) - (V))(V(0) - (V))}: its time 
integral / °° C(t')dt' gives an exact measure of the diffusion 
coefficient. This time integral can also be estimated by as- 
suming an exponential decay C(t) = T r /Mexp(— t/r): in 
this case D T = T r r/M. The asymptotic decay, in fact, 
is always exponential, since the process is Markovian by 
construction. This is checked in the insets of Fig. 5. 
From Eqs. (11)-(13) it can be deduced that, in the large 
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e 2 

Fig. 5: The diffusion coefficient as a function of e 2 calculated in 
two different ways (see text)for Q2 = 0.3 and oli = 0.9. Inset: 
a) The self-diffusion d 2 (t) = ((d(t) - (d(t))) 2 } of the piston ver- 
sus time; b) The self-correlation function of the piston-velocity 
versus time. The dashed line corresponds to the theoretical 
prediction of D T (see text) for e 2 = O.Of , a? = 0.9 and r c = 2.5. 



M/m limit the dynamics of the piston is analogous to an 
Ornstcin-Uhlenbcck process with an effective friction 1/r 
and noise intensity 2T r /r, with r = r c /[e 2 (3 + a^)] and 
t c is the average collision time. From such an expression 
it is seen that the diffusion coefficient does not depend on 
the mass and very little on the inelasticity. For smaller 
values of M/m, on the contrary, the diffusion coefficient 
increases, as displayed in Fig. 5. A growing discrepancy 
between D T and D is observed in this limit, due to the fact 
that the first part of the correlation decay is not exactly 
exponential. 

Conclusion. — We have discussed a simple example 
of noise rectification when two identically driven granular 
gases are separated by a moving piston. The difference 
of piston inelasticites between the two faces, induces a 
stationary drift and the appearance of non-Gaussian fluc- 
tuations, which become highly asymmetric when inertia 
is reduced. A description of the dynamics in terms of 
the first three moments of the velocity distribution is suf- 
ficient to predict most of the physics of this system. It 
is quite rare to reach such a detailed knowledge of a sys- 
tem with asymmetric non-Gaussian fluctuations: this is 
achieved, here, thanks to the empirical observation that 
the third cumulant of the distribution is dominant over 
cumulants of higher order, a fact that is not guaranteed 
in general [17]. We are also confident that the assumption 
used, mainly that the two gases on the sides of the pis- 
ton are very dilute and driven at high frequency, can be 
reproduced in the laboratory and could be at the base of 
ratcheting mechanisms observed in real granular materi- 
als. 
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A.l. — The coefficients of the Eqs. (8)-(10) can be written as 
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where = (1 + a 2 )/2 and A(7j) = y/r)/(l + rje 2 ). 
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